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Abstract. The shell model Monte Carlo (SMMC) method is a powerful technique for 
calculating the statistical and collective properties of nuclei in the presence of correlations 
in model spaces that are many orders of magnitude larger than those that can be treated 
by conventional diagonalization methods. We review recent advances in the development and 
application of SMMC to mid-mass and heavy nuclei. 


1. Introduction 

Most microscopic studies of heavier nuclei are based on mean-field methods such as density 
functional theory. However, important correlation effects beyond the mean field can be missed. 
The configuration-interaction (Cl) shell model approach accounts for correlations and shell 
effects but conventional diagonalization methods are limited by the size of the many-particle 
model space. The auxiliary-field Monte Carlo (AFMC) method enables calculation of thermal 
observables and ground state properties in model spaces that are many orders of magnitude 
larger than model spaces that can be treated by conventional diagonalization methods. In 
nuclear physics this method is also known as the shell model Monte Carlo (SMMC) method U,!2]. 
SMMC has proven to be a powerful method in calculating statistical and collective properties 
of nuclei 1312 0 017]. 

Here we review recent advances in the development and application of SMMC to nuclei. In 
particular, we describe a method to circumvent the sign problem that arises from the projection 
on an odd number of particles and its application to calculate pairing gaps in mid-mass nuclei 
and level densities of odd-mass isotopes. We also discuss the application of SMMC to describe 
the microscopic emergence of collectivity in heavy rare-earth nuclei and present results for the 
collective enhancement factors of level densities in such nuclei. 

2. Auxiliary-field Monte Carlo (AFMC) method 

The AFMC method utilizes the Hubbard-Stratonovich transformation [8], in which the 
imaginary-time propagator e~^ H of a nucleus that is described by a Hamiltonian H at inverse 


temperature /3 = 1/T is represented by a functional integral of one-body propagators of non¬ 
interacting nucleons over external auxiliary fields 
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Here G a is a Gaussian weight and U a is the propagator of non-interacting nucleons moving in 
time-dependent auxiliary fields <t(t). The thermal expectation value of an observable O is 

Tr(e~P H O) _ f D[a]W a $ a (0) a 
Tr e~P H J D[a]W„<lG 



where W a = G a \Tr U a \ is a positive-definite weight function and = Tr U a /\ Tr U a \ 
is the Monte Carlo sign function. The quantity ( 0) a = Tr (OU a )/TicU a is the thermal 
expectation value of the observable for a given configuration a of the auxiliary fields. In 
SMMC, we choose samples cr k distributed according to W a and estimate the observables from 
(0)t ~ J2k(0)a k ^a k / ■ AFMC often suffers from a sign problem at low temperatures. 

However, there are classes of good-sign interactions that are free of the sign problem and are 
composed of the dominant collective components of effective nuclear interactions [9j. Small 
bad-sign components can be treated by the extrapolation method of Ref. [2]. 


3. Circumventing the odd-particle sign problem 

Applications of SMMC to odd-even and odd-odd nuclei have been hampered by a sign problem 
that originates from the projection on a odd number of particles. Such a sign problem occurs at 
low temperatures even for good-sign interactions and makes it difficult to estimate accurately 
the ground-state energy of an odd-particle system. A breakthrough was a method we introduced 
to calculate accurately the ground-state energy of an odd-particle system that circumvents this 
sign problem m- It is based on the scalar imaginary-time single-particle Green’s functions 

G u{r) = X^ Ta ^ m ( r ) a tm( 0 )> ■ ( 3 ) 

m 

Here T denotes time ordering, u = ( nlj ) denotes the radial quantum number n, orbital angular 
momentum l and total spin j of the single-particle levels, and m is the magnetic quantum 
number of j. In an asymptotic regime of sufficiently large r but t <C /3, the Green’s functions 
of an even-even nucleus A have the form G u (r) ~ e P\- E A A ± l ) £ g . s .(A)]M ; where Ej(A± 1) is the 
lowest energy of spin J = j in the neighboring even-odd nuclei A ± 1 and -E g . s . is the ground-state 
energy of the nucleus A. Thus Ej(A ± 1 ) can be determined from the slope of the calculated 
lnG^ and the ground-state energy of the even-even nucleus. The ground-state energy of the 
even-odd nucleus is then found by minimizing Ej(A ± 1 ) over j . 

We used this method to calculate neutron pairing gaps in mid-mass nuclei from the second- 
order energy difference in neutron number A n = (— ) N [E(Z, N +1) + E(Z, N — 1) — 2 E(Z, N)]/ 2, 
where E(Z, N ) is the ground-state energy of the nucleus with Z protons and N neutrons. Fig. [I] 
shows neutron pairing gaps for mid-mass isotopes. The calculations were carried out in the 
complete fpgg /2 shell with the interaction of Ref. [3]. The SMMC results (solid circles) are in 
overall good agreement with the experimental values (open circles). 


4. Level densities 

SMMC has been successful in calculating nuclear state densities p(E x ) versus excitation energy 
E x in the presence of correlations. The density calculated in SMMC is the state density, in 
which each level with spin J is counted 2J + 1 times. However, often the density measured in 
the experiments is the level density, in which each level with spin J is counted once. 





Figure 1. Neutron pairing gaps A„ as a function of mass number A in families of fpg g/ 2 -shell 
isotopes. The gaps calculated in SMMC (solid circles) are compared with the experimental gaps 
(open circles). Adapted from Ref. [TO] . 


We recently introduced a method to calculate the level density directly in SMMC El- 
Denoting by pm the level density at given value M of the angular momentum component J z , 
the level density p is given by p = pm =o for even-mass nuclei and by p = Pm= 1/2 for odd-mass 
nuclei. J z projection can be carried out exactly in SMMC as discussed in Ref. [5]. 

Level densities were recently measured by proton evaporation spectra in a family of nickel 
isotopes 59 ~ 64 Ni [T2] . We calculated these level densities [13] within the fpgg /2 shell using the 
Hamiltonian of Ref. j]3]. The calculation of the ground-state energies of the odd-mass nickel 
isotopes, and hence their excitation energies, was made possible by the method discussed in 
Sec. [ 3 J Fig. [|] shows the level densities of 59-64 Ni. The SMMC results are in close agreement 
with the level densities determined from various experimental data sets. 



Figure 2. Level densities p in a family of nickel isotopes. The SMMC level densities (solid 
circles) are compared with level densities determined from various experimental data sets: proton 
evaporation spectra (open squares and quasi-continuous lines), level counting data (histograms) 
and neutron resonance data when available (triangles). Taken from Ref. [13] . 


5. Microscopic emergence of collectivity in heavy nuclei 

Heavy nuclei are known to exhibit various types of collectivity. Typically, nuclei near shell 
closure are spherical and display vibrational collectivity, while mid-shell nuclei are strongly 










deformed and display rotational collectivity. Different types of collectivity are well described by 
empirical models but much less is understood microscopically. Here we use SMMC to describe 
the emergence of collectivity in heavy nuclei within the framework of the Cl shell model approach. 


5.1. Thermal signatures of collectivity 

A given type of collectivity is usually identified by its associated spectrum. While SMMC is 
a powerful method for calculating thermal and ground-state properties in very large model 
spaces, it does not provide detailed spectroscopic information. To overcome this difficulty, we 
have identified a thermal observable whose low-temperature behavior is sensitive to the type of 
collectivity. Such an observable is (J 2 )t> where J is the total angular momentum. Assuming an 
even-even nucleus with a rotational or a vibrational ground-state band and an excitation energy 
E 2 + of the first 2 + state, we have [6J [7] 
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Figure 3. Left: (J 2 )t hr 162 Dy. The SMMC results (solid circles) are compared with fits 
to the rotational model (solid line) and vibrational model (dashed-dotted line). Adapted from 
Ref. [6]. Right: (J 2 )t hr 148 Sm. The SMMC results (open circles) are compared with fits to the 
vibrational model (solid line) and rotational model (dashed line). 


We carried out calculations for rare-earth nuclei using the 50 — 82 shell plus I/ 7/2 for protons, 
and the 82 — 126 shell plus 0h 11 / 2 , 159/2 f° r neutrons. The Hamiltonians we used are given in 
Refs. [6, 7]. The left panel of Fig. [ 3 ] shows (J 2 )t for the deformed rare-earth nucleus 162 Dy. 
The SMMC results are compared with fits to both the rotational and vibrational models. 
The agreement of the SMMC results with the rotational model confirms that we are able to 
describe the rotational character of lb2 Dy in the framework of a truncated spherical shell model 
approach. The fit to the rotational formula gives E 2 + = 84.5 ± 8.9 keV, in agreement with the 
experimental result of E 2 + = 80.6 keV. The right panel of Fig. [ 3 ] shows (J 2 )t for the spherical 
rare-earth nucleus 148 Sm. Here the SMMC results follow closely the vibrational formula fit with 
E 2 + = 538 ± 31 keV, in agreement with the experimental value of E 2 + = 550 keV. 

5.2. The crossover from vibrational to rotational collectivity 

Heavy nuclei are known to exhibit a crossover from vibrational to rotational character as 
the number of neutrons increases from shell closure towards mid-shell. In a mean-field 
approximation, this is described by a phase transition from spherical to deformed nuclei. We 
have studied two such families of rare-earth isotopes, 148 ~ 154 Sm and and 144 ~ 152 Nd [7J. 










Figure 4. (J 2 )t as a function of temperature T for a family of even neodymium isotopes. The 
solid lines describe the experimental values extracted using Eq. ([ 5 ]). For 150 ~ 152 Nd we only use 
the discrete levels in Eq. ([ 5 ]) (dashed lines). Adapted from Ref. [7J. 


Fig. [ 4 ] shows (J 2 )t for a family of even neodymium isotopes. The open circles are the SMMC 
results. We observe a crossover from a soft response to temperature in 144 - l4fi Nd. typical of a 
vibrational nucleus, to a rigid response in 152 Nd, typical of a rotational nucleus. The solid lines 
are extracted from experimental results using 

(J 2 )t = + + dE x p(E x ){J 2 ) Ex e- E */ T ^ , (5) 

where Z{T) = + l)W Bi / T + J^ N dE x p(E x )e~ Ex ^ T is the corresponding experimental 

partition function. Here Ei are the measured low-lying energy levels with spin Jj, which we 
assume to be complete up to to certain threshold energy En- Above the threshold energy, we 
use a backshifted Bethe formula for the level density p(E x ) whose parameters are determined 
from a fit to level counting data at low energy and neutron resonance data at the neutron 
resonance threshold. We use in Eq. ([ 5 ]) the spin cutoff model estimate (J 2 )e x = 3er 2 , where 
cr 2 = IT /ft 2 is calculated assuming a rigid-body moment of inertia I. We observe an overall 
good agreement between the SMMC values of (J 2 )r and the experimentally determined values. 

6. Collective enhancement factors of level densities in heavy nuclei 

SMMC state densities in rare-earth nuclei were found to be in good agreement with experimental 
data |14| . Collective states in these nuclei lead to enhancement of their state density, which is 
described by a collective enhancement factor K. The decay of K with excitation energy is one of 
the least understood topics in the modeling of level densities m, and it is often parameterized 
by empirical formulas. We define a collective enhancement factor by K = psmmc/VhfBi where 
/Osmmc is the total SMMC state density and phfb is the level density calculated in the finite- 
temperature Hartree-Fock-Bogoliubov (HFB) approximation. The latter describes the density 
of intrinsic states and thus the above ratio is a measure of the enhancement in the number of 
states due to collective states that are built on top of the intrinsic states. 

Fig.0 shows the calculated enhancement factor K as a function of excitation energy E x for 
the even neodymium isotopes 144 ~ 152 Nd. We also show by arrows the energies of the various 
thermal phase transitions that occur in the HFB approximation. The thin arrows correspond to 
the pairing transitions (for both protons and neutrons), while the thick arrows describe the shape 
transitions. The spherical nuclei 144_146 Nd do not support rotational bands, so the enhancement 
factor K is due to vibrational states alone. We observe that this vibrational enhancement factor 
decays to ~ 1 at excitation energies above the proton and neutron pairing transitions. The 
heavier neodymium isotopes 148_152 Nd are deformed in their ground state. We observe that K 
has a local minimum above the pairing transitions, which we attribute to the interplay between 









Figure 5. Collective enhancement factor K (open circles) of the even neodymium isotopes 
144 ' 152 Nd versus excitation energy E x (see text for the definition of K ). The thin arrows 
correspond to the proton and neutron pairing transitions, while the thick arrows describe the 
shape transitions. Adapted from Ref. m- 


the decay of vibrational collectivity and the enhancement due to rotational collectivity. The 
plateau that follows the local minimum describes the enhancement from rotational states. This 
rotational enhancement decays in the vicinity of the shape phase transition. 


7. Conclusion 

The SMMC is a powerful method for calculating statistical and collective properties of nuclei in 
the framework of the Cl shell model in very large model spaces. We reviewed a recent method to 
circumvent the odd-particle sign problem and its application to the calculation of pairing gaps 
and level densities of odd-mass nuclei. We discussed recent applications of SMMC to heavy 
rare-earth nuclei and in particular the microscopic emergence of collectivity in such nuclei. 
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